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The hydrodynamic (hydro) model has been extensively applied to heavy ion data from the rela- 
tivistic heavy ion collider (RHIC). Results are interpreted to conclude that a dense QCD medium 
nearly opaque to most partons, a strongly-coupled quark-gluon plasma (sQGP), is formed in 
more-central Au-Au collisions. The sQGP may have a very small viscosity ("perfect liquid"). 
However, other analysis methods provide contradictory evidence. Two-component analysis of 
single-particle hadron spectra reveals a spectrum hard component consistent with a parton frag- 
ment distribution described by pQCD which can masquerade as "radial flow" in some hydro- 
motivated analysis. Minimum-bias angular correlations reveal that a large number of back-to- 
back jets from initial-state scattered partons with energies as low as 3 GeV survive as "minijet" 
hadron correlations even in central Au-Au collisions, suggesting near transparency to partons. 
In this talk I present methods by which almost all spectrum and correlation structure, even in the 
most-central Au-Au collisions at 200 GeV, can be described quantitatively by pQCD calculations. 
The evolution of nuclear collisions is apparently dominated by parton scattering and fragmenta- 
tion even in the most-central A-A collisions, albeit the fragmentation process is strongly modified. 
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1. Introduction 



The intention of this work is to test the extent to which perturbative QCD (pQCD) can describe 
more-central A-A collisions at RHIC. Is a hydrodynamic (hydro) description necessary, or even 
allowed by data? Detailed arguments are provided in Refs. [jl|, ^, |3| |]], with related material on 
hydro interpretations of azimuth quadrupole structure in Refs. []|, ||, ^, g]. 

I begin with the two-component spectrum model for p-p collisions. Then I review the phe- 
nomenology of fragmentation functions (FFs) from LER HERA and Fermilab. I describe calcula- 
tions on that basis of perturbative QCD (pQCD) fragment distributions (FDs). What was described 
this morning as "fragmentation functions" I distinguish from true FFs which are distributions con- 
ditional on parton energy. Fragment distributions can be calculated by folding FFs with a parton 
spectrum. I introduce a parton "energy-loss" model [^] to provide FD calculations which can de- 
scribe measured fragmentation evolution with A-A collision centrality. The last part of the talk 
extrapolates beyond the single-particle system to describe a method for converting jet angular cor- 
relations into fragment yields and spectra. It is then possible to determine quantitatively the minijet 
contribution to the A-A final state. What emerges is a comprehensive pQCD description of RHIC 
nuclear collisions up to central Au-Au. 



2. Visualizing fragmentation: Conceptual consequences of plotting formats 

I first consider the impact of plotting-format choices on physical interpretations of data: For in- 
stance, to what extent does a given plotting format favor hydro over fragmentation interpretations? 
Figure [T] (first panel) shows a conventional plotting format for fragmentation functions, beautiful 
LEP data from OPAL at 91 GeV that Yuri showed you. In this format the featured structure at larger 
x p actually represents a small fraction of the fragment yield — less than 10% of the fragments — that 
which can be described by DGLAP evolution. Some people are quite interested in those details. 
The structure at upper left in the first panel (small x p ) is typically ignored at RHIC. 

Figure [I] (second panel) shows the same data replotted on normalized rapidity u: rapidity 
variable y = ln{(E + p)/m 7[ } divided by the equivalent measure y max = \n{Q/m n ) for the parton. 
We find that FF data so plotted can be described by a beta distribution (solid curve) to the error 



limits of the data down to zero momentum fllPp. The solid curve is not a theoretical description, it 
is a phenomenological description of measured FFs. DGLAP applies toward the right and MLLA 
applies near the peak, but the beta distribution accurately describes all fragment data from parton 
momentum down to zero momentum (in both panels). 

Figure [l] (third panel) shows conventional spectrum ratio Raa for five centralities of 200 GeV 
Au-Au collisions plotted on transverse rapidity y t = ln{(m t + p t )/m n } ^ and p-p data from 



Ref. [|1 1|. Raa is defined as the ratio of an A-A p t spectrum to a reference p-p spectrum divided 
by the (Glauber) number of binary collisions. The p-p (N-N) reference p#jv is in this case the two- 
component model described in the next section. The conventional emphasis is on p t above 6 GeV/c 
(upper scale) which is 4.5 on y t (lower scale). Information about fragmentation below 6 GeV/c is 
strongly suppressed by Raa, presenting a misleading picture. 

Alternatively, one can extract from the same spectra "hard components" Haa, which are most 
relevant to fragmentation, and replot them as ratio yaa 0]- In Fig, [I] (fourth panel) yaa reveals 
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Figure 1: First: Fragmentation function for 91 GeV e + -e~ collisions Jl2| ] plotted on conventional variable 
x p . Second: The same data plotted on normalized rapidity u. The solid curve is a beta distribution. Third: 
Conventional spectrum ratio for five centralities of 200 GeV Au-Au collisions (bold curves) and p-p data 
(solid dots). The thin solid curves are linear superposition references. Fourth: The same data plotted in the 
form of hard-component ratio taa, revealing large enhancements at smaller p t corresponding to suppression 
(jet quenching) at larger p t . 



for the first time that hard-component variations at 0.5 GeV/c (large enhancement) are exactly 
complementary to variations at 10 GeV/c (suppression). Centrality evolution in the two places is 
strongly correlated. We should not claim an understanding of fragmentation (e.g. jet quenching) 
until the entire fragmentation picture is acknowledged. 



3. Two-component model of spectra and correlations 

The two-component model has several manifestations at RHIC. I refer to the two-component 
spectrum model as first developed in Ref. JTl|]. The basic physical model is similar to that in 
PYTHIA, but the details are determined by data phenomenology. In p-p collisions the "soft compo- 
nent" refers to longitudinal fragmentation of projectile nucleons by soft-Pomeron exchange leading 
to diffractive dissociation. The "hard component" refers to large-angle scattered parton fragmen- 
tation from a minimum-bias parton spectrum, possibly by hard-Pomeron exchange. The fragment 
hadron distribution extends in principle down to zero momentum. 

Figure |2| (first panel) shows p t spectra from non-single-diffractive (NSD) p-p collisions at 



= 200 GeV [11]. The spectra correspond to ten observed (uncorrected) multiplicities n c f, in 
one unit of tj. Corrected multiplicity n c h ~ 1n c h- The spectra have been normalized by "soft" 
multiplicity n s determined iteratively by a limit process. The data are plotted on transverse rapidity 
y t with pion mass representing unidentified hadrons. The spectra are then functions of y t and n c h. 

Evolution with h c h is simple. The spectra can be represented accurately by a Taylor expansion 
in fi c h which has only two terms. The "constant" term So is a fixed function of y t , and the "linear 
coefficient" Hq is also a fixed function of y t , both independent of n c f,. The two-component model 
for p-p collisions with soft and hard multiplicities n s + «/, = n c h is then 

1 ldnchi * ch) -S (y t ) + n -^H (y t ), (3.1) 



n s {n ch )y t dy t n s (h ch ) 

Factor tih/n s is observed to vary as ah c h. So(y t ) is by definition the limiting spectrum as n c f, — > 
and has the form of a Levy distribution on m t . By subtracting So(y t ) from each of the spectra in the 
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Figure 2: First: Spectra for ten multiplicity classes of 200 GeV NSD p-p collisions normalized to soft- 
component multiplicity n s [p"T|]. So(y t ) is the "soft component" limiting case for n c /, — > 0. Second: Spectrum 
hard components inferred from the same spectra by subtracting the soft component from all spectra. Third: 
Two-component model of the p-p spectra. Ho(y t ) is the hard-component model function. Fourth: Average 
hard component for NSD p-p collisions (data points) with pQCD "fragment distribution" (solid curve). 



first panel we obtain the residuals in the second panel. The form is independent of multiplicity and 
well described by the solid curves representing fixed form Ho(y t ), a Gaussian plus QCD power- 



law tail on transverse rapidity y t Qll|]. The third panel shows the model in Eq. (|3JJ) which can 
be compared with data in the first panel. For comparisons with A-A spectra (below) we define 
S pp = (l/yt) dn s I dy t with reference model n s Sq and similarly for H pp tit, Ho- 

In Fig. ^| (fourth panel) the points (spectrum hard component) represent an average of the 
p-p hard components in the second panel, each scaled to the multiplicity density corresponding to 
NSD p-p collisions. The solid curve represents a calculated pQCD fragment distribution discussed 
below. 

The corresponding two-component model for per-participant-pair A-A spectra is 

2 1 dn c h , . , 
-r" = Smbt) + v H AA(y t ,v) (3.2) 

tlpartyt dy t 

= S NN (y t ) + VrAA(yt,v)H NN (y t ), 



where Snn (~ S pp ) is the soft component and Haa is the A-A hard component (with reference 
Hjw ~ Hpp) ^ [fy]. Ratio yaa = Haa/H^n is an alternative ratio measure to nuclear modification 
factor Raa- Centrality measure v = Inuruay / ^ n participant estimates the mean projectile-nucleon path 
length in A-A collisions. We are interested in the evolution of hard component Haa or ratio yaa 
with A-A centrality. For the A-A two-component model the spectrum soft component remains by 
hypothesis unchanged and scales as the number of participant pairs n part /2. For Glauber linear 
superposition of p-p (N-N) collisions (GLS reference) spectrum hard component Haa — > H^(y t ) 
would also remain unchanged modulo the factor v relative to participant scaling. In real A-A col- 
lisions HAA(yt,b) changes relative to GLS reference fl^(y t ), representing "medium modification" 
of parton fragmentation. 

Study of HAA(yt,b) reveals evolution of fragmentation with centrality pi]. pQCD can be 
used to calculate equivalent fragment distributions which can be compared directly with measured 
Haa l@l ■ FD calculations require a combination of measured fragmentation functions, a pQCD- 
predicted parton spectrum and a theoretical model of FF modification in A-A collisions. Several 
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sections below describe methods to calculate FDs and present interpretations of spectrum hard- 
component evolution based on theory-data comparisons. 

4. Fragmentation functions 

We require a phenomenological representation of measured fragmentation functions. As noted, 
FFs are commonly plotted on momentum fraction x p = p fragment / P jet or = m (!/•*/>)■ Alterna- 
tively, one can define rapidity y = ln{ (E+p)/ m K } with pion mass adopted for unidentified hadrons, 
where p is the fragment total momentum appropriate for most FF data from e + -e~ collisions. 
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Figure 3: First: Measured fragmentation functions (points) from e + -e~ collisions for three CM (dijet) 
energies JT^ , [l3| ] plotted on rapidity y with /3 -distribution parametrizations (solid curves). Second: The 
same FF data and curves scaled to unit-normal distributions and plotted on normalized rapidity u. Third: 
Parametrization of the ensemble of e + -e~ FFs (surface) over a large energy range. Fourth: FF ensemble 
parameterization for p-p (p-p) collisions. 

Figure || (first panel) shows FFs from e + -e~ collisions at three energies plotted as D(x, Q 2 ) o 
D(y : y max ) = 2dn c h/dy on y [ 12, 13]. To good approximation the FFs are self-similar: both widths 
and amplitudes scale linearly with y max = \rv{Qjm^) = ln(2Ej et /m n ) JTOj]. The jet fragment mul- 
tiplicity n c h j(y max) is then approximately proportional to y 2 max . Deviations from that trend arise 
mainly from the running of a s . Because the FFs are nearly self-similar we can convert them to 
a universal form by renormalizing both the FF amplitude (to unit-normal) and the rapidity [to 
u = (y — y min) / (ymax — Vmm)> where y m j n ~ 0.3]. We then obtain the data in the second panel 
represented by a beta distribution (solid curves) to the error limits of the data: D(y,y max ) — > 
1-n c hj(ymax)fi(u;p,q) (where p,q are parton-energy-dependent parameters). Only the best FF data 
reveal any deviations from perfect y max (Q 2 ) scaling. The self similarity is the dominant aspect of 
DGLAP evolution. The simplicity is not apparent unless FFs are plotted on y. If the universal beta 
distribution is transformed back to individual jet energies we obtain the solid curves in the first 
panel and the surface in the third panel, which accurately describes all FFs above 3 GeV partem 
energy (2 = 6 GeV) and down to zero hadron fragment momentum [10]. Each vertical slice of the 
surface plot is the FF for a particular parton energy. 

CDF p-p FF data (not shown) exhibit significant differences from e + -e~ FFs. The p-p FF 
ensemble is represented in the fourth panel as modified e + -e~ FFs. Part of the difference is due to 
an imposed cone radius which should exclude some low-momentum fragments. However, the p-p 
data suggest a real reduction relative to e + -e~ FFs. The evolution with energy scale of p-p FFs is 
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also anomalous: there is a saturation of the FF amplitudes at larger jet energies compared to LEP 
FFs. FF universality may not be a valid assumption given the energy trend of those data [gj. 



5. Parton spectrum 

Next we require a pQCD parton spectrum. In Fig. || (first panel) the solid curve is a power-law 
spectrum with energy cutoff inferred by working backward from 200 GeV p-p and A-A spectrum 
hard-component data according to procedures described below. The bold dotted curve is an ab 
initio pQCD calculation [M]. The two spectra agree quantitatively near 3 GeV (y max ~ 3.8) where 
almost all scattered partons appear. The good agreement, arising from two independent approaches 
to parton spectrum determination, is significant. Parton spectrum details at larger parton energies 
are less important for minimum-bias spectrum and correlation analysis. 




E. d (0.eV) 10 j 




Figure 4: First: Parton spectra inferred from this analysis for p-p collisions (solid curve) and central Au- 
Au collisions (dash-dotted curve) compared to an ab-initio pQCD theory result (bold dotted curve [|l4|l). 
Second: Fragment distribution (solid curve) compared to p-p hard-component data (points). Dotted curves 
correspond to ±10% change in parton spectrum cutoff energy about 3 GeV. Third: Curves in the first panel 
plotted on a linear scale to illustrate that almost all partons (gluons) appear near 3 GeV. Fourth: Comparison 
of parton spectrum inferred in Ref. [||] (solid curve) with ab initio pQCD calculations (bold dashed and 
dot-dashed curves) and event-wise reconstructed jet spectra (data points) [Il5 16 1. 



Given the power-law approximation to the parton spectrum it is not apparent from pQCD 
where the effective spectrum cutoff should be. Figure |] (second panel) shows calculated FDs from 
the procedure described below with spectrum cutoffs at 3 ± 0.3 GeV (solid and dotted curves) 
compared to the p-p spectrum hard-component data introduced above. The comparison establishes 
an empirical 3 GeV spectrum cutoff. Figure ^ (third panel) shows the power-law spectrum (with 
cutoff) on a linear scale, demonstrating that almost all scattered partons appear at the cutoff energy. 

Saturation-scale (SS) arguments support a cutoff at 1 GeV (light dotted curve in the first 
panel) fll!]]. Given the approximate power-law dependence the difference in the total parton yield 
for the two cutoffs is a factor 30-50 in the initial parton (mainly gluon) density. The 1 GeV SS 
cutoff is based on an argument derived from initial-state parton densities in nucleons. Considered 
as a quantum-mechanical process parton scattering and fragmentation to charged hadrons depends 
not only on the initial-state parton density but also on the final-state hadron density of states. If 
there is no final state for a given parton scattering the transition is not allowed. The effective cutoff 
should then depend on the available density of hadronic final states at a given parton energy scale. 
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Figure || (fourth panel) shows data from a UA1 analysis of energy clusters in EM calorimeter 
data (solid dots) leading to inference of "minijets," with spectrum extending down to 5 GeV (later 
amended to 3-4 GeV after background subtraction) [|l^]. More recent data from the STAR collab- 
oration (open circles) extending down to 5 GeV (4 GeV without background contribution) are also 
shown flip]]. The parton spectrum inferred phenomenologically from RHIC p-p data (solid curve) 
includes a cutoff near 3 GeV which is consistent with the UA1 observations of 1985 and with the 
STAR reconstructed-jet spectrum. The undershoot of the solid curve at larger parton energy may be 
due to oversimplified modeling of p-p FFs. If the p-p FF saturation mentioned above is included the 
inferred parton spectrum should be even closer to pQCD theory and jet spectra. The dash-dotted 
curve is a spectrum from Ref. [17]. 



6. pQCD folding integral and fragment distributions 



The pQCD folding (convolution) integral used to calculate fragment distributions is 
d 2 n h e(AT})/2 f°° , n , ^dGdi 



dydrj cr NSD ATj 4jr Jq 



dy m a X D{y,y n 



iijet 



dy„ 



(6.1) 



where D{y,y max ) is the FF ensemble for a specific collision system (e + -e~, p-p, A-A, in-medium 
or in-vacuum), and dGdijet / dy max is the parton spectrum [||]. The perturbative object is the parton 
spectrum; the nonperturbative object is the measured FF ensemble. The folding integral then pro- 
duces a prediction for the observed hadron spectrum hard component. d 2 ni l /dydr\ is the predicted 
hadron FD from parton pairs scattered into angle acceptance Atj. Efficiency factor e 6 [1,2] in- 
cludes the possibility that the second jet of a dijet also falls within Atj . At^ ~ 5 is the effective 471 
tj interval for scattered partons [[jj], and Onsd ~ 36 mb is the cross section for NSD p-p collisions, 
both for v^wv = 200 GeV. 





Figure 5: First: pQCD folding-integral argument for e + -e~ FFs. Second: e + -e~ FD (solid curve) and 
p-p hard-component reference from Ref. [ pi] ] (dash-dotted curve). Third: Folding-integral argument for p-p 
FFs. Fourth: p-p FD (solid curve), p-p hard-component data (solid dots) and reference (dash-dotted curve). 

Figure || (first panel) shows the folding-integral argument D(y,y l1iax )dadijet/dy max for e + -e~ 
FFs and 3 GeV parton spectrum cutoff. The fragment distribution is then the projection of the 
2D histogram onto fragment rapidity y. Figure |] (third panel) shows the argument for FFs from 
p-p collisions with cutoff on fragment rapidity higher than for e + -e~ collisions. Figure |5] (second 
panel) shows the FD projection (solid curve) compared to a model of the measured spectrum hard 
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component in NSD p-p collisions (dash-dotted curve). The FD from e + -e~ FFs is adopted below 
as the reference for all measured hard-component distributions. In general, fragment distributions 
(FDs) from theory are compared to spectrum hard components (HCs) from data. Figure ||] (fourth 
panel) shows the FD (solid curve) for p-p collisions compared to the HC (points) from NSD p-p 
collisions [|l(i|]. That comparison established the 3 GeV partem spectrum cutoff [Q]. 

According to these calculations the most significant variations in fragmentation and the largest 
fragment yields appear below 2 GeV/c, which reveals & fundamental logical problem in the con- 
ventional RHIC approach to data analysis and interpretation. The spectrum interval below 2 GeV/c 
(y t = 3.3), described as the "soft" region of the hadron spectrum, is conventionally assigned to 
hydro models. The interval above 6 GeV/c is conventionally assigned to "hard processes" (par- 
ton scattering and fragmentation) described by pQCD. The source of the logical difficulty lies in 
confusing theoretical limitations on pQCD descriptions of fragmentation functions with the theo- 
retical ability to describe fragment distributions in terms of measured FFs. It is the parton spectrum 
which must be described perturbatively, not fragmentation functions. By imposing an unjustified 
constraint on pQCD descriptions of fragmentation the great majority of hadron fragments is sur- 
rendered to hydro interpretations. 

7. Parton "energy loss" and medium-modified fragment distributions 

We next require a model for medium modification of fragmentation functions. Figure ^ (first 
panel) shows "medium modified" fragmentation functions (Borghini-Wiedemann, BW) achieved 
by altering certain splitting functions in the parton cascade [Q]. The solid and dashed curves (vac) 
are parametrizations of e + -e~ FFs based on the beta distribution which describe FF data within 
their uncertainties down to zero fragment momentum [[k]]. The dotted and dash-dotted curves 
(med) are FFs modified to match the BW prescription simply by changing the parameter q in the 
beta distribution which, by construction, conserves the parton energy. There is suppression of 
larger-momentum fragments and consequent enhancement of smaller-momentum fragments p]. 




y = ln[(E + p)/ng y max y y 



Figure 6: First: e + -e~ FFs for two energies unmodified [ |10| | (solid and dashed curves) and modified to 
emulate parton "energy loss" []|] (dash-dotted and dotted curves). Second: Modified e + -e~ FF ensemble. 
Third: Medium-modified FD from e + -e~ FFs (solid curve) compared to in-vacuum FD (dotted curve). 
Fourth: Medium-modified FD from p-p FFs (solid curve) compared to in-vacuum FD (dotted curve). 

The modified e + -e~ FF ensemble is shown in the second panel. The locus of modes (white 
dashed curve) is shifted to smaller fragment rapidities y (compare with Fig. |3| third panel). There 
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is a similar result for the p-p (p-p) FF ensemble. By inserting modified FFs into the pQCD folding 



integral Eq. ( p.l| ) we obtain modified fragment distributions (solid curves) compared to unmodified 
FDs (dotted curves) in the third and fourth panels At larger p t there is "jet suppression" 
(parameter q is adjusted to match spectrum data there). At smaller p t there is corresponding "jet 
enhancement," a new aspect of the fragmentation problem. The low-p f enhancement is large for 
e + -e~ FFs but negligible for p-p FFs. That difference becomes important in A-A collisions. 



8. Evolution of fragmentation with centrality 

We can now make a direct comparison between calculated pQCD FDs and measured spectrum 
HCs. Figure ^ (first panel) shows measured pion spectra from 200 GeV Au-Au for five centralities 
plotted on pion rapidity (dark solid curves) [fxlj . Snn is the common soft component inferred as 
the limiting spectrum for centrality measure v — > 0. The solid dots are the hadron spectrum from 
NSD p-p collisions [11]. The dash-dotted curve is hard component H^n inferred from the n c h 



dependence of the p-p spectrum. 

Figure [7] (second panel) shows Haa extracted from the Au-Au spectra (bold solid curves) 



according to Eq. ( p.2[ ): we subtract the same soft component from spectra for five centralities 
and divide by v. Although the soft component certainly dominates spectra below 1 GeV/c the 
systematic uncertainty in Haa is manageable at least down to 0.5 GeV/c (y t = 2). The fixed soft 
component plus the inferred hard components describe the original spectrum data exactly. Relative 
to GLS reference H^n (bold dashed curve) there is suppression at larger p t and enhancement at 
smaller p t in more-central Au-Au collisions. Detailed systematic study of HC evolution with Au- 
Au centrality reveals that with increasing centrality i) p-p FFs transition to e + -e~ FFs, ii) FFs 
become "medium modified" and iii) there is a 50% increase in the dijet cross section due to a 10% 
reduction in the effective parton spectrum cutoff energy (3 GeV — > 2.7 GeV). The bold dotted 
curves in the second panel show pQCD FDs calculated according to those systematic trends [0]. 



1 p,(GeV/c) 10 , 
200 GeV Au-Au 
v < 2.5 




Figure 7: First: Pion spectra (dark solid curves) for five centralities of 200 GeV Au-Au collisions. Second: 
Hard-component centrality evolution in Au-Au collisions [Jl|]. Large enhancements at smaller y t accompany 
suppression at larger y t . Third: Hard-component (HC) ratios relative to an ee-vacuum reference for Au-Au 
collisions below the sharp transition. Fourth: HC ratios above the sharp transition revealing major changes 
in HC structure. 

To obtain a more differential view of HC evolution we define ratios r xt = H^/Href, where 
xx denotes the collision system and H re f is not inferred from p-p collisions. Instead, H re f — > 
H ee -vac is defined by an FD constructed according to Eq. (6.1) using the parton spectrum inferred 



9 



Partem fragmentation in nuclear collisions 



Thomas A. Trainor 



from p-p collisions [Q] and in- vacuum FFs from e + -e~ collisions JTQ|] - Figure [7| (third panel) shows 
r xx data for the HC from p-p collisions (solid dots) and from more-peripheral Au-Au collisions 
(bold solid curve). Reference curve r^N (dashed) is obtained from the FD for N-N collisions using 
p-p FFs (~ //aw)- Dash-dotted reference curve r ee is obtained from the FD for in-medium modified 
e + -e~ FFs corresponding to central Au-Au collisions (parton spectrum cutoff reduced by 10%). 

Figure [7] (fourth panel) shows r„ data (bold solid curves) for HCs from more-central Au-Au 
collisions. We observe a sharp transition in HC evolution at v = 2.5, with qualitatively different 
behavior below and above the transition. Below the transition we observe strong suppression of 
the HC at smaller p t relative to what is expected for in-vacuum e + -e~ FFs (third panel). Above 
the transition, in more-central Au-Au collisions, we observe a strong enhancement at smaller p t 
complementing the suppression at larger p t observed with conventional ratio Raa- The HC central- 
ity dependence at 0.5 GeV/c corresponds in detail to that at 10 GeV/c. The number of particles 
at smaller p t resulting from parton fragmentation is much greater than that at larger p t (consistent 
with approximate parton energy conservation). The copious low-;?, hadron fragments should be 
accommodated in any theoretical description of A-A collisions. 



9. Fragment yields from jet angular correlations 



We can infer fragment yields from /^-integrated jet angular correlations via factorization of 
the measured jet-correlated pair density to obtain the per-jet fragment multiplicity. We integrate 
Eq. ( |6.1[ ) over fragment rapidity y on both sides to obtain dn^ / df] , the per-participant-pair fragment 
density on T], in terms of jet frequency f(b) and mean jet fragment multiplicity n c hj{b) 

f, d 2 n h f e(Ari)o dijet (b) \ f 1 [°° da dijet \ 

/ dy = < > < 7TT- / dy max n ch j (y max ,b) — } (9.1) 

J dydri [ o NSD kr\ An J {o di j et {b) J dy max J 

= f(b)n chJ (b). 

n c hj(b), averaged over the minimum-bias parton spectrum, is effectively the mean fragment multi- 
plicity for partons near the parton spectrum cutoff (~3 GeV). We can infer n c hj(b) from jet angular 
correlations for A-A centrality b. Jet frequency f(b) = {\/n},i n )dnj(b)/dr] is the number of jets 
per unit tj per NSD N-N collision estimated from pQCD. The argument of Gdij e t(b) admits the 
possibility that the N-N dijet cross section may depend on A-A centrality [||]. 



9.1 Jet angular correlations 



2D angular autocorrelations on difference variables t]a = T]i — T]2 an d <Pa = (j>\ — 02 evaluated 



near mid-rapidity retain all angular correlation information [|18[]. 2D correlations can be constructed 
for the p t -integral minimum-bias case or with specific p t cuts on one or both particles in a pair. 2D 



angular correlation histograms are formed for p-p collisions [ ]19j , [20J and several (typically 11) 
centrality classes of A-A collisions [|2l|]. Figure [8] (left panels) shows 2D histograms for peripheral 
(v = 1.4, ~ p-p) and mid-central (v ~ 4.8) 200 GeV Au-Au collisions. The correlation structure 
has three main features: a same-side (SS, 0a < ^/2) 2D peak at the origin, an away-side (AS, 
<Pa> 7l/2) ridge approximately uniform on T]a and described by dipole cos (0a — 7r) in more-central 



A-A collisions, and non-jet azimuth quadrupole cos(20a)- Angular correlations in Ref. [21] are 
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reported in the per-particle form Ap / y/p re f [18]. The SS jet peak is modeled by a 2D Gaussian 



Apss 



= p (b)f(r} A ,(j) A ,b) =A 2D exp (-t?1/2ctJ) exp (-0a/2<tJ) , 



(9.2) 



where j' 2 is a sibling/mixed pair ratio and po(b) is the single-particle 2D angular density at mid- 
rapidity. 

In Fig. |8] (left panels) the fitted non-jet quadrupole has been subtracted. The SS peak can 
be interpreted as intrajet correlations and should include all hadron fragment pairs from jets that 
survive partonic and hadronic rescattering. The AS ridge can be interpreted as internet correlations 
from back-to-back scattered partons. The most probable p t for minimum-bias jet-correlated pairs 
is 1 GeV/c in p-p collisions, consistent with the mode of the spectrum hard component. 




Figure 8: First: 2D angular correlations from more -peripheral (v ~ 1.4) 200 GeV Au-Au collisions. 
Second: 2D angular correlations from more-central (v ~ 4.8) 200 GeV Au-Au collisions. Third: Same-side 
2D peak amplitude (solid curve) compared to GLS reference (A-A transparency, dashed curve). Fourth: 
Same-side 2D peak widths (solid curves) compared to GLS references (dashed curves). 



SS 2D peak properties vary strongly with centrality. Figure M (right panels) shows SS peak pa- 



rameter variations with centrality parameter v Q2l|]. We observe a sharp transition in SS peak prop- 
erties at v ~ 2.5 in 200 GeV Au-Au collisions which corresponds to the spectrum hard-component 
systematics noted in Sec. ||. Below the transition the peak properties follow Glauber linear superpo- 
sition of p-p structure as expected for transparent A-A collisions. Above the transition the SS peak 
amplitude increases rapidly relative to GLS, there is strong elongation on tj and slight narrowing 
on (j) . The narrowing on azimuth is inconsistent with parton energy-loss models based on multiple 
scattering. It is notable that even in more-central Au-Au collisions (e.g. Fig. ||, second panel) 
the SS peak is well-described by a 2D Gaussian, and the AS peak on azimuth is an undistorted 
dipole SS peak systematics appear to correspond to jets, but what do those correlation trends 
imply for single-particle yields in the final state? To answer that question we convert (factorize) 
two-particle jet correlations to obtain the equivalent in single-particle hadron fragment yields and 
spectrum hard components. 

The SS peak volume is by hypothesis the number of jets in the angular acceptance times the 
number of fragment pairs per jet, which allows us to factorize the SS jet peak. To convert from 
jet angular correlations to parton fragment yields and spectra requires four steps: i) angle-average 
the SS 2D peak on (t]a, <Pa) over the 4D angular acceptance on (t]i ,772,^1,^2) to obtain mean pair 
ratio j 2 (b), ii) estimate the mean pQCD jet number per event nj(b) within the 77 acceptance, iii) 
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calculate the mean per-jet fragment multiplicity n c hj{b), iv) combine those elements to infer jet 
fragment yields/spectra as HC 2D densities on (rj,<p). The result of step i), the average of the SS 



2D peak described by Eq. ( |9.2| ) over the angular acceptance, is shown in Fig. |9| (first panel) in the 
form po(b)j 2 (b). The dashed curve shows the result for a 471 acceptance. 



9.2 Jet properties from jet correlations 

Fig. ^(second panel) shows apQCD estimate of jet frequency f(b) (solid curve) on A-A cen- 
trality measure v. The datum marked by the open symbol was inferred from a two-component 
analysis of spectra for 200 GeV p-p collisions with At] = 1 [[TT|]. The increase near v = 2.5 cor- 
responds to the observation that the parton spectrum cutoff energy drops by about 10% near the 
sharp transition, leading to an approximate 50% increase in the dijet cross section [0]. 
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Figure 9: First: Same-side 2D peak averaged over the angular acceptance (solid curve) and corresponding 
result for 4k acceptance (dashed curve). Second: pQCD estimate of jet frequency f{b) in Au-Au collisions 
(solid curve). Third: Corresponding mean jet number per event vs Au-Au centrality (solid curve) and corre- 
sponding binary collision scaling (dashed curve). Fourth: Mean per-jet fragment multiplicity inferred from 
trends in previous panels within the angular acceptance (lower curves) and in An (dashed curve). 

Figure || (third panel) shows the corresponding number of jets rij (b) = nj,j n Ar) f(b) for Arj = 2 
within the STAR TPC. The sharp transition in jet correlation structure and spectra occurs at about 
the point in Au-Au collision centrality (v ~ 2.5, upper hatched region) where the number of jets 
becomes significantly greater than one. The dijet cross section is observed to increase by about 
50% at the same centrality. The dashed curve (GLS) is N-N binary collision scaling. 

We now combine estimated jet frequencies with measured jet angular correlations to infer 
mean jet fragment multiplicities. Figure ^ (fourth panel) shows the mean multiplicity derived from 
jet angular correlations. The assumption that the number of jet-correlated pairs is equal to the 
number of jets times the mean number of fragment pairs is expressed in the first line of 



nj(b)n 2 chJ (b)=n 2 ch (b)j 2 (b) 



(9.3) 



n c h,j(b) = n ch (b) y / ' j 2 (b)/tij(b). 



The second line expresses the required factorization, where n c h(b) = 2nAr\ po(b) (charged-particle 
multiplicity in the angular acceptance) and j 2 (b) are measured quantities. Mean per-jet fragment 
multiplicity n c h j is thus inferred from correlation data and a pQCD jet number hypothesis. The 
fragment multiplicity for untriggered jets (mainly 3 GeV minijets) is something between 2 and 
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4 for p-p collisions, increasing to about 6 in central Au-Au, within the angular acceptance. The 
dashed curve is what the fragment multiplicity would be with 4% acceptance. In more-central Au- 
Au collisions jets are elongated on tj, and part of the jet falls outside the TPC angular acceptance. 
Ironically, n^ j is more uncertain in p-p than in central Au-Au collisions. The several curves 
and hatched regions illustrate the systematic uncertainty in the multiplicity estimate. We can now 
calculate the fragment density in Eq. ( JO| ) and the minimum-bias jet contribution to the final state. 



9.3 Minijet (minimum-bias parton fragment) contribution to the final state 



Figure 10 (first panel) shows spectrum hard component HaaQj) (solid curve) inferred from jet 



angular correlations according to Eq. ( |9.1| ) 



(9.4) 



The open point is an estimate from Ref. Qllp . The solid points are derived from the "total hadrons" 
data in Fig. 15 (left panel) of Ref. [jlj|. Multiplying through by v/2n gives the first line of 



VHAA(b) 



2 ( ,,n chJ {b) 
-rijib) 



L part 



L part 



27TATJ 



(9.5) 



The second line incorporates the second line of Eq. (£3j) and the definition of single-particle den- 
sity po(b). vHAA{b) is the hard component in the two-component spectrum model of Eq. (fO[). 
Figure [l0| (second panel) shows the two-component particle yield Snn + vHaa (b) predicted by 
measured jet angular correlations (bold solid curve). Soft component snn is by hypothesis fixed at 
~ 0.4 [2D density on (17, 0)] for all A-A centralities. The solid points are the "total hadrons" data 
in Fig. 15 (left panel) of Ref. M\ divided by 2n. 
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Figure 10: First: Inferred p t -integrated hard-component yield vs centrality. The open symbol is an estimate 
from Ref. pTTj], Second: Total charged-particle yield vs centrality estimated from two-particle jet correlations 
(solid curve) and a two-component representation of measured single-particle data (dash-dotted line). Third: 
Single-particle 2D angular density po vs centrality. Participant scaling is indicated by the dashed curve. 



The dash-dotted curve is the Kharzeev-Nardi (KN) approximation to the per-participant 2D 
density (2/n part )po(b) measured in more-central Au-Au collisions p2[]. Figure [10| (fourth panel) 
shows charged-hadron density po(b), the 2D angular density on (tj,0) assumed for this analysis. 
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The solid curve is the KN model po(b) = (n part /2)p pp {\ +x{v — 1)} with p pp = 0.4 and x = 0.09 
for Au-Au at 200 GeV p2fl. The KN description matches minimum-bias data in more-central 
collisions but fails in more-peripheral collisions where corrected yield data are sparse. For more- 
peripheral collisions we expect a GLS trend extrapolated from p-p collisions (dotted line in second 
panel). The sharp transition in jet angular correlations near v = 2.5 explains the deviation. 

Equation ( |9"3| ) (second line) implies that y/nj(b) j 2 (b) is the fractional hadron yield from par- 
ton fragmentation (minijet fractional yield). The first factor in the radicand is obtained from pQCD 
(relative systematic uncertainty < 50%). The second factor is from measured jet angular corre- 
lations (relative uncertainty small). Combining angular correlation measurements and a pQCD 
estimate of jet number we find that one third of the hadronic final state in central Au-Au collisions 
at 200 GeV is associated with resolved jet correlations (relative uncertainty < 20%). 

10. Summary 

Hydro-motivated analysis of RHIC data tends to interpret the large hadron fragment contribu- 
tion below 2 GeV/c in terms of flow phenomena. The role of parton scattering and fragmentation 
in nuclear collisions is minimized. Its pQCD description is artificially restricted to small regions of 
phase space. In contrast, model-independent analysis of spectrum and correlation structure reveals 
new fragmentation features quantitatively described by pQCD over the full momentum range. 

Hard components extracted from p, spectra are now identified as single-particle manifesta- 
tions of minimum-bias parton fragmentation in nuclear collisions. Spectrum hard components 
correspond quantitatively to minimum-bias jet angular correlations (minijets). pQCD fragment 
distributions calculated by folding a minimum-bias parton energy spectrum with a parametrization 
of measured fragmentation functions accurately describe the measured hard components. 

Modification of fragmentation functions in more-central A-A collisions can be modeled by 
adjusting a single parameter in the FF parametrization (beta distribution) consistent with rescaling 
specific QCD splitting functions. The reference for all fragmentation in nuclear collisions is the 
FD derived from measured in-vacuum e + -e~ FFs and the parton spectrum for p-p collisions. 

Relative to the reference the spectrum hard component for p-p and peripheral Au-Au collisions 
is found to be strongly suppressed for smaller fragment momenta. At a specific point on centrality 
the Au-Au spectrum hard component transitions to strong enhancement at smaller momentum and 
suppression at larger momentum, described by FDs derived from medium-modified e + -e~ FFs. 

Minimum-bias jet (minijet) correlations have been converted to absolute fragment yields which 
are found to comprise approximately one third of the final state in central 200 GeV Au-Au colli- 
sions. Those results indicate that almost all large-angle scattered partons down to 3 GeV parton 
energy survive as jet manifestations in the final state, albeit with some modification. 

Novel effects in A-A collisions may be related to strong color fields established among scat- 
tered energetic partons, an elaboration of three-jet events in LEP collisions. Hadron fragment 
structure may directly reflect the large-scale color field geometry as a manifestation of local parton- 
hadron duality (LPHD). Newly-interpreted spectrum and correlation systematics, correctly associ- 
ated with parton fragmentation, suggest evolution of the color-field geometry in nuclear collisions. 
We conclude that pQCD calculations should be applied to all aspects of spectrum and correlation 
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data in order to discover what is truly novel in RHIC collisions. We find that perturbative QCD can 
describe a large part of RHIC collision evolution - and hydro interpretations are questionable. 
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